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ABSTRACT 

We have observed with the 30-m IRAM telescope, the CSO telescope and the ISO 4 satellite (Infrared 
Space Observatory) the rotational lines of CO at millimeter, submillimeter and far infrared wavelengths 
in the direction of C-rich stellar objects at different stages of evolution : CRL 2688 (a very young 
Proto-Planetary Nebula), CRL 618 (a Proto-Planetary Nebula), and NGC 7027 (a young Planetary 
Nebula). Several changes in the longwave emission of CO and other molecules are discussed here in 
relation with the degree of evolution of the objects. In the early stages, represented by CRL 2688, the 
longwave emission is dominated by CO lines. In the intermediate stage, CRL 618, very fast outflows 
are present which, together with the strong UV field from the central star, dissociate CO. The released 
atomic oxygen is seen via its atomic lines, and allows the formation of new O-bearing species, such as 
H2O and OH. The abundance of HNC is enhanced with respect to HCN as a result of the chemical 
processes occurring in the photo-dissociation region (PDR). At this stage, CO lines and [01] lines are 
the dominant coolants, while the cooling effect of [CI J] is rising. At the Planetary Nebula stage, NGC 
7027, large parts of the old CO AGB material have been reprocessed. The spectrum is then dominated 
by atomic and ionic lines. New species such as CH + appear. Water has probably been reprocessed in 
OH. 

Subject headings: infrared: stars — line: identifications — planetary nebulae: individual (CRL 2688, 
CRL618, NGC7027) — stars: abundances — stars: carbon — stars: evolution 



1. INTRODUCTION 

Solar type stars remain in the asymptotic giant branch 
(AGB) for about ~ 10 6 — 10 7 years, losing mass through 
stellar winds (M ~ 10~ 4 — 10 -7 M Q /yr, ejection velocity 
~ 5 — 25 kms -1 , see, e.g., Loup et al. 1993). During this 
phase the ejected material progressively forms an expand- 
ing circumstellar envelope (CSE). Molecular species are 
easily formed in the innermost regions of the CSE. Some 
of these molecules aggregate onto dust grains which un- 
der the action of radiation pressure accelerate and push 
the remaining gas producing an expanding dusty molecu- 
lar envelope. UV photons from the interstellar radiation 
field dissociate gas-phase molecules in the external lay- 
ers of the envelopes allowing new chemical reactions and 
the production of new molecular species. These chemi- 
cal changes have been very well studied through observa- 
tions and models in the case of IRC+10216, the prototyp- 
ical AGB carbon-rich star (Glassgold 1996; Cernicharo, 
Guelin & Kahane 2000 and references therein). Addi- 
tional changes in the composition of the CSE do occur 
when the central object starts its evolution towards the 
white dwarf stage. A large UV field arises from the cen- 
tral (much hotter) star and at the same time, high velocity 
winds appear and interact with the AGB remnant (Cer- 
nicharo et al. 1989; Kwok 2000). Objects in this phase are 
called proto-planetary nebulae (hereafter PPNe). During 
this period (1000-2000 years; Bujarrabal et al. 2001), the 
almost spherical symmetry typical of the AGB stage disap- 



pears and is replaced by a complicated geometrical struc- 
ture (most PPNe are elliptical, bipolar or quadrupolar, 
Zuckerman & Aller 1986; Frank et al. 1993). The chemi- 
cal composition is also strongly affected by these processes: 
O-bearing molecules can be formed in C-rich objects (Her- 
pin & Cernicharo 2000), and complex organic molecules, 
which are not observed in IRC+10216, are efficiently pro- 
duced (Cernicharo et al. 2001 a&b). 

In order to better understand this evolution, we present 
in this paper a comparative study of the millimeter, sub- 
millimeter and far-IR CO line emission from 3 objects rep- 
resenting different stages of this fast transition: CRL 2688, 
a very young PPN, CRL 618, a PPN, and NGC 7027, 
a young planetary nebula (PN). The different excitation 
coSnditions of the observed CO lines allow to probe differ- 
ent layers in the CSEs. In particular, we study the CSE 
remnant wind and the higher velocity winds. The CSE 
component is the "normal" AGB wind, and the wind com- 
ponent refers to a higher velocity wind which is probably 
not present during the AGB phase and which contributes 
to changing the morphology from spherical to bi-polar or 
other non-spherical geometries. The observations are pre- 
sented in section 2. The representativity of each selected 
source is discussed in section 3. Section 4 is devoted to 
the analysis of the spectra of IRC+10216 as a reference of 
the AGB phase. In section 5 we discuss the wind struc- 
ture of the young PPN CRL 2688. Section 6 is devoted to 
CRL 618 and section 7 to the analysis of NGC 7027 ob- 
servations. An overall comparative discussion of the four 
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objects is then presented in section 8. 

2. OBSERVATIONS AND DATA ANALYSIS 

The IRAM-30m observations of the CO 1-0 and 2-1 
lines were performed in September 2000 using four SIS re- 
ceivers with high (>0.90) image sideband rejection. The 
pointing and focus were checked using either Saturn or 
Jupiter, using the continuum emission of the sources and 
was found to be accurate to within 2". The typical system 
temperatures during these observations were 250 K (CO 
1-0) and 400 K (CO 2-1). The observations used the 
wobbler system which provides very flat baselines. The 
spectrometers were a filtcrbank (at 115 GHz) and an au- 
tocorrelator (at 230 GHz) set to cover a band of 512 MHz 
with frequency resolutions of 1 and 1.25 MHz respectively 
(velocity resolutions of 2.60 and 1.63 kms -1 respectively). 

The CO 3—2, 4—3, 6—5 and 7—6 observations were per- 
formed with the 10.4 m telescope of the Caltech Submil- 
limeter Observatory at the summit of Mauna Kea (Hawaii) 
in September and October 2000 using three different re- 
ceivers based on helium-cooled SIS mixers operating in 
double-sideband mode (DSB). The typical system temper- 
atures during these observations were 750 K (CO 3—2) and 
5000-6000 K (CO 4-3, 6-5 and 7-6). For further details 
about the receivers see Kooi et al. (2000) and references 
therein. The backends were two different acousto-optic 
spectrometers. One of them covers a bandwidth of 500 
MHz with 1024 channels, and the other one covers 1.5 
GHz with 2048 channels thus providing respectively the 
following velocity resolutions (in kms -1 ): 146. 38/^ and 
219.57/^, with v in GHz. The instantaneous band of the 
receivers, however, is ~ 0.9 GHz, thus limiting the validity 
of the data from the 1.5 GHz AOS to that range. 

Since some of the CSO observations were performed 
shortly after a new 790—920 GHz receiver was installed 
for the first time at the telescope, its pointing and focus 
were carefully determined during previous runs earlier in 
the year. Then, the pointing of this and the other re- 
ceivers was checked during our observations using Jupiter 
and Saturn first and then the own CO emission of each ob- 
ject. Taking into account that this emission is somewhat 
extended, the pointing is considered accurate to within 
4" (CO 3-2, 4-3 lines) and 2" (CO 6-5 and 7-6 lines). 
The observations were also performed by wobbling the sec- 
ondary reflector by 60". Proper corrections were finally 
applied to the data to compensate for the use of only one 
load (ambient temperature) in the chopper-wheel standard 
calibration of the CSO that relies upon the Penzias-Burrus 
(1973) approximation. 

For the comparative discussion that we make further in 
this paper, fits of the CO mm/submm observations (shown 
in Fig. 1) of the three objects have been performed to de- 
rive the flow velocity of the different components and the 
integrated intensity ratio Ri between the emission of the 
CSE remnant and that of the wind for each transition. 
As angular resolution changes when observing at differ- 
ent frequencies and using different telescopes, the relative 
contribution to the bulk of the detected emission from the 
CSE remnant (extended) and the high velocity winds (un- 
resolved) will be different. Hence, Ri outlines the wind 

3 The ISO Spectral Analysis Package (ISAP) is a joint development 
institutes are CESR, IAS, IPAC, MPE, RAL and SRON. 



excitation of the different flows across the envelope (cf. 
Tables 1,2 and 3). 

The ISO LWS Fabry-Perot data (A/AA - 9500) of CRL 
2688 were obtained during orbits 366 and 369. The in- 
tegration times are 3400 and 5800 sec. The spectral res- 
olution is 0.015 ^m. The LWS grating mode (43-196.7 
/im, A/AA ~ 200) has a spectral resolution of 0.29 ^m 
for the 43—93 /im range and 0.6 /jm for the 80—196 /im 
range. LWS01 grating spectra for CRL 2688 and CRL618 
were taken during orbits 21 and 688 respectively (total on- 
source times of 1218 sec and 5360 sec, respectively). The 
J=18-17 line of HCN was observed towards CRL 2688 
during orbit 154 with the LWS grating mode and an in- 
tegration time of 914 sec (see Fig. 4). Finally, the LWS 
grating spectra of NGC7027 taken by ISO during orbits 21, 
342, 349, 356, 363, 377, 537, 552, 559, 566, 579, 587, 594, 
601, 706, 713, 720, 727, 734, 741, 755, 762, 769, 776 and 
783 were averaged. The total on-source time was 53409 
sec. All data were processed following pipeline number 9 
and treated with ISAP 3 to remove glitches and fringes. 



3. REPRESENTATIVELY OF THE SOURCE SAMPLE 

In the next sections we try to identify and discuss the 
changes that ocurr through the transition of a C-rich AGB 
star to the PN stage. The aim of this work is to trace 
the evolution of the circumstellar chemistry in the late- 
AGB and post- AGB phases. This study is based on the 
analysis of the spectral differences observed in 4 objects: 
IRC+10216, CRL 2688, CRL 618 and NGC 7027. We 
thus have to relay upon the assumption that IRC+10216 
will evolve into an object similar to the Egg Nebula, and 
then to an object similar to CRL 618 which will finally 
evolve to a planetary nebula that looks like NGC 7027. 
The ideal, but unfeasible approach, would be to observe 
a single object evolving through the post- AGB phase. In 
addition, our sample objects should be "typical" in the 
sense that each one should represents the averaged char- 
acteristics that the bulk of objects in the same stage of 
evolution. 

However, CRL 2688, CRL 618 and NGC 7027 seem to 
be quite unusual objects in some properties, extreme cases 
among all carbon-rich AGB and post-AGB objects. All 
arc very bright in molecular lines, but while IRC+10216 is 
a relatively close source (110-150 pc, Groenewegen et al. 
1998, Crosas & Menten 1997), the other three objects are 
5-10 times further away (Acker et al. 1992, Bujarrabal et 
al. 1994, Sahai et al. 1998). This may twist the compari- 
son between IRC+10216 and the other sources, as distance 
influences line fluxes and profiles. 

Recent works estimate a < 2 M Q progenitor mass for 
IRC+10216 and place its core mass to 0.6 M Q (Kahanc 
et al. 2000). Theoretical calculations suggest a progeni- 
tor main sequence mass of ~ 3 M Q for both CRL 2688 
and CRL 618 (Speck, Meixncr and Knapp 2000), and 3-4 
Mq for NGC 7027 (Sahai et al. 2001). IRC+10216 there- 
fore, seems likely to evolve to the PN phase rather more 
slowly than the others, and this may affect the chemistry of 
the envelope. Anyway, several parallelisms exist between 
these objects: dense, optically thick dusty circumstellar 

the LWS and SWS Instrument Teams and Data Centers. Contributing 
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Fig. 1. — Spectra for the 12 CO J=l-0, 2-1, 3-2, 4-3, 6-5 and 7-6 transitions for CRL 2688, CRL 618 (no J=4-3 observation) and NGC 
7027. The main beam temperatures are in K and the velocity in kms -1 . Baselines were applied. The 1—0 and 2—1 lines are IRAM 30m 
observations; the others were made with the CSO; lines are smoothed. For each line, the upper part of the caption shows the whole profile, 
and the lower part shows the high velocity emission with an expanded vertical scale. Note that the emission present on the red part (at 141 
kins -1 ) of the CO (3—2) emission from CRL 618 is the HC3N (J=38— 37) line emission (at 345.609 GHz). The conversion factors Jy/K are 
respectively 4.3, 4.2 and 40 for the 1—0, 2—1 observations and rest of observations. 
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envelopes and shell-like density enhancements observed as 
incomplete arcs (Men'shchikov et al. 2001, Mauron & Hug- 
gins 1999). These circumstellar concentric arcs, extremely 
pronounced in IRC+10216, are also detected in NGC 7027 
(see Kwok, Su & Stoesz 2001). More generally, it has been 
proposed that further evolution of CRL 2688 will lead to 
morphologies similar to that observed in NGC 7027 (spe- 
cially the ionized nebula and H2 region; Cox et al. 1997). 

The Egg Nebula is quite anomalous compared to more 
usual carbon-rich post-AGB objects. Its infrared spec- 
trum is very different from that of any other carbon-rich 
post-AGB object, since it appears to be optically thick 
even at 30 /im (Omont et al. 1995). More common are the 
most extreme dusty carbon-stars such as the small group 
of sources discussed by Volk et al. (2000). 

CRL 618 is quite similar to CRL 2688 in that both have 
massive molecular envelopes and central bipolar outflows 
(Meixner, Fong & Sutton 2001). However, CRL 618 shows 
an unusual rich circumstellar chemistry that could be re- 
lated to its advanced stage of evolution (Herpin & Cer- 
nicharo 2000). 

NGC 7027 is quite unusual among planetary nebulae in 
that it has a relatively massive circumstellar envelope be- 
yong the ionized region. Its PDR is one of the densest and 
warmest in the PN sample studied by Liu et al. (2001). 
In fact, the ionized region is known to be relatively low 
in mass (< 0.2 solar masses, smaller than the typical val- 
ues of 0.3-0.5 Mq of the ionized region of PN) while the 
surrounding mostly molecular gas totals 2-4 solar masses. 

In conclusion, the objects selected here are rather pe- 
culiar in some properties that may influence their post- 
AGB evolution. However, a comparative study of their 
molecular envelopes could be useful to identify the dis- 
tinct chemistries induced by the general physical changes 
produced by the different degrees of evolution, and this is 
the aim of this paper. 

4. IRC+10216 : THE PROTOTYPICAL AGB CARBON-RICH 

OBJECT 

To introduce this study, it is useful to present the main 
features of a typical object at the end of the AGB stage. 
IRC+10216 is the brightest AGB carbon rich star. Its esti- 
mated distance is 110—150 pc (Groenewegen 1997, Crosas 
& Menten 1997). It is believed to be in the final stage of 
red giant evolution and exhibits a bipolar reflection neb- 
ula (Skinner, Meixner & Bobrowsky 1998). The temper- 
ature of the central object is 2000 K, and the expansion 
velocity of the CSE is 14.5 kms -1 (Skinner, Meixner & 
Bobrowsky 1998; Cernicharo et al. 2000). IRC+10216 has 
an extended CSE where more than 50 molecular species 
have been detected (see line surveys by Cernicharo et al. 
1996, Cernicharo et al. 2000). This object has a rich 
carbon chemistry and many of the species detected in its 
molecular envelope are carbon chain radicals which are 
formed in the external layers of the CSE (Cernicharo et 
al. 1996; Cernicharo et al. 2000). The innermost regions 
of the envelope are dominated by a chemistry at thcrmody- 
namical equilibrium (Tsuji 1973), where CO, HCN, C2H2, 
SiO are formed. As shown by Cernicharo et al. (1996), 
the far-IR spectrum consists of strong dust emission and 
several molecular emission lines from 12 CO (J=14— 13 to 
J=39-38), 13 CO ([ 12 CO / 13 CO]~ 45), HCN (J=18-17 



to J=48-47), H 13 CN and vibrationally excited HCN. The 
HCN lines have intensities similar to those of CO, and 
[HCN/CO]= 1/10, proving that HCN is the main coolant 
in the innermost region of the CSE. Only one ion, HCO + , 
has been detected (Lucas & Guelin 1999) and there is no 
[CI I] line emission. The mid- and near-infrared spec- 
tra arc dominated by the ro-vibrational lines of C2H2 
and HCN (Cernicharo et al., 1999). The most abun- 
dant molecules through the CSE are CO, HCN and C2H2 
(Fuente, Cernicharo & Omont 1998, Cernicharo et al. 
2000). 

The molecular observations of this object are well ex- 
plained by a photochemical model (Glassgold 1996) con- 
sidering that, due to the low effective temperature of the 
central star, the only UV photons available in the outer 
CSE come from the interstellar radiation field. 

5. A YOUNG PPN: CRL 2688 

CRL 2688 (the Egg Nebula) is a very young PPN (it 
has left the AGB phase probably only about 100 years 
ago) with an effective temperature around 6600 K (spec- 
tral type F5, Justtanont et al. 1997) and luminosity close 
to 1.8 10 4 Lq (Jura & Kroto 1990). The non-detection 
of [OI] and [CII] (Cox et al. 1996) indicates that radia- 
tion from the central star has not yet formed a PDR. The 
molecular gas is in expanding fragmented shell structures, 
and shocks are thought to heat it (Cox et al. 1996). The 
AGB CSE expansion velocity is 20 kms -1 . CO emission 
has been seen as far as 42" from the central object. The 
weaker 13 CO emission extends over 20" (Yamamura et al. 
1996). Sahai et al. (1998) detected a moderate velocity 
wind (hereafter MVW) of 40 kms -1 with an angular ex- 
tension of 10". High velocity gas emerging in two different 
directions from the central star has been detected by Cox 
et al. (2000) in the CO 3=2—1 emission line. The size 
of the emitting region is 4" with extensions in both the 
north-south and east-west directions. Cox et al. (2000) 
have shown that this high velocity gas is closely related to 
the H2 emission at 2 /ito. 

5.1. Low- excitation CO lines 

The profiles of all low excitation CO lines (see Fig.l) 
are characteristic of optically thick emission. However, 
while the 3=1— 0, 2—1, and 3—2 lines probably trace the 
extended AGB remnant envelope of CRL 2688, higher J 
transitions are more sensitive to hotter gas from the inner 
layers of the envelope and hence of lower spatial extent. 
As J increases, the line opacities in the AGB CSE decrease 
and the emerging profiles change accordingly. 

From the CO observations of CRL 2688 and the line 
parameters given in Table 1 we can distinguish two main 
outflows. The main one, the AGB CSE remnant, is cen- 
tered at —35.4 kms -1 and covers around 20 kms -1 (in 
good agreement with Young et al. 1992). The second out- 
flow has a moderate velocity extent, ~50 kms -1 , and is 
centered at —60 kms" 1 . The intensity of the red wing of 
this MVW is weaker than its blue counterpart in the CO 
lines. We believe that this behavior is produced by AGB 
CSE remnant absorption of the MVW red wing. 

An absorption dip is detected in the blue part of the 
mm and submm CO lines at ~— 55 kms -1 , just at the 
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Fig. 2. — LWS/FP observations of 12 CO lines in CRL 2688 (left panel) and average spectra of all observed lines (right 
panel). The intensity scale corresponds to the continuum normalized flux rescaled to the values obtained with the LWS 
gratings for the same lines. The wavelength of each transition is indicated at the top left corner of each panel. The right 
panel shows the average spectra of the CO lines. Dashed lines show the fit produced by our LVG model as described in 
the text (with a center velocity precision of ±15 kms -1 ). 



terminal velocity of the AGB CSE. It suggests that the 
low-velocity wind absorbs radiation from the inner and 
faster wind. Young ct al. (1992) found that this absorp- 
tion has a velocity width of ~1 kms -1 in the low- J lines. 
We found the same value for the J=6— 5 and 7—6 CO lines. 
This would suggest that the microturbulence through the 
AGB CSE is nearly constant and equal to 1 kms -1 . 

The ratio Rj goes from 0.08 (J=l-0) to 0.24 (J=7-6), 
stressing the importance of the MVW emission in the 
intcrmediatc-J CO lines. The dilution of the MVW emit- 
ting region in the telescope beam will change with J. In 
addition, line opacities for the AGB CSE remnant and the 
MVW will also change differently with J. Hence, above 
a given rotational transition, the emission from the AGB 
CSE remnant will be weaker than that from the MVW. 
Observations of high excitation lines are thus needed in 
order to follow up the CO emission from the inner enve- 
lope. 

5.2. Far-IR emission 

Far-IR observations reveal the inner and hotter regions 
of the CSE where important chemical and physical pro- 
cesses are changing the structure of the AGB CSE. Rota- 
tional lines of 12 CO (J=14-13 to 26-25) are clearly de- 
tected while the emission from 13 CO is very weak. HCN 
is clearly detected in the AOT LWS02 grating spectrum 
of its J=18— 17 transition (see Fig. 4, line detected with 



S/N ~ 4; all the sub-peaks are smoothed noise). Note that 
some CO lines present broad wings at the grating resolu- 
tion (see J=19-18, J=20-19, J=21-20) that are probably 
instrumental artefacts resulting from the defringing of the 
LWS data. 

Several high-J CO (from J=16— 15 to 26—25) lines were 
observed in Fabry-Perot mode with the LWS (see Fig. 2 for 
the J=16-15, J=19-18, J=24-23 and J=26-25 lines). 
The fluxes were divided by the continuum and rescaled to 
the values obtained for the same lines with the LWS grat- 
ing. To derive more valuable information about kinematics 
and because the excitation and opacity is rather similar in 
all these lines, we computed an average spectrum (better 
S/N ratio) of all the individual lines (see right panel in Fig. 
2). Study of this average CO LWS/FP spectrum reveals 
a much broader line than what we could expect from the 
AGB CSE remnant and centered at —65 kms -1 . It could 
prompt that all the emission is coming from the MVW. 
In order to estimate the CO emission at high-J from the 
AGB CSE remnant and from the MVW, we have fitted 
two gaussians to the average spectrum of Figure 2. We 
obtain velocities of — 70±10 kms -1 and — 35±10 kms -1 
(see on the same figure the modelled emission from our 
Large Velocity Gradient (LVG) calculations -discussed be- 
low) . These velocities agree, within the FP absolute wave- 
length calibration, with the MVW and AGB CSE veloci- 
ties respectively. Obviously, higher spectral resolution ob- 
servations are required to confirm both components (i.e., 
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Fig. 3. — Continuum subtracted LWS spectra of CRL 2688 (bottom caption), CRL 618 (middle caption) and NGC 7027 (top caption). The 
result of our models is shown by the continuous grey line. The lines of CO , CO , HCN, H2O , CH+ and OH arc indicated by arrows while 
those of HNC in CRL 618 arc indicated by vertical lines (from J=22-21 at 150.627 fan to J=17-16 at 194.759 urn). The [CI I] and [0/7/] 
lines are not included in our models; the plots indicate gaussian fits to these features. 
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HERSCHEL observations with HIFI). In the following we 
assume that both components emit in the high-J lines of 
CO and we will try to derive the relative contribution of 
each one. 

Obviously, the average spectrum does not directly pro- 
vide information on the column density of CO. In order 
to derive this parameter we have fitted the CO J=16— 15 
emission observed with the FP (the line with the best S/N 
ratio) with the two velocity components derived from the 
average spectrum. We then extended the result to the 
other lines (Fig. 2, left panel). Globally, the lines are 
well fitted by our model, with column densities of 85% 
(feature at —70 kms -1 ) and 45% (feature at —35 kms -1 ) 
with respect to the values retrieved from fitting the grating 
spectrum. 

The red counterpart of the emission at —70 kms -1 is 
not detected. Our lines suggest that one part of the emis- 
sion (at —70 kms -1 ) comes from a wind with an outflow 
velocity of about 35 kms -1 , probably the MVW of Sahai 
et al. (1998). As the —70 kms -1 emission comes from a 
region where large shocks occur (interaction with the AGB 
remnant envelope), and the column density is thought to 
be larger, this emission appears stronger. 

The other component (at the stellar velocity) is consis- 
tent with the wind velocity derived from our CO mm and 
submm observations (see Table 1). The widths of the emis- 
sions in each CO line show that the velocity dispersion is 
probably important in the excitation region. 

In order to model the observed emission we have used 
an LVG code which divides the emitting source into sev- 
eral components with different solid angles (defined by 8). 
For each component we have fitted the column density of 
the molecule under consideration, the H 2 volume density 
and the temperature. Hence, the molecular abundance 
will depend on the layer thickness which cannot be de- 
rived in this retrieval scheme. While these assumptions 
could be correct for the high velocity wings and the PDR 
(CRL 618, NGC 7027), the component corresponding to 
the AGB remnant is poorly approximated as the whole 
CSE is reduced to a layer with a given average column 
density, temperature and density. This AGB remnants has 
been modelled in detail by Bujarrabal & Alcolea (1991). 

The same basic physical parameters have been applied 
to 13 CO and HCN and we then only varied the column 
densities (see parameters in Table 4). For all the objects, 
model components were selected according to previous ob- 
servational and theoretical works and our millimeter and 
submillimeter observations. The particular temperatures, 
solid angles, etc. are also derived from these studies, even 
if in some cases we adapt the parameters to correctly fit 
the data. Of course, column densities for each layer are 
sensitive to the excitation temperature. 

According to other authors and our study of the low-J 
CO line emission, the adopted geometry model and phys- 
ical parameters of each component are : 

(i) An inner region representing the fast wind from the 
star (Sahai et al. 1998) and modelled as a layer of 8= 0.3" 
, expansion velocity of 100—200 kms -1 , temperature of 
800 K, and a density of 5 10 7 cm" 3 . The 12 CO column 
density (N) is 2 10 18 cm" 2 and the 13 CO and HCN abun- 
dances relative to 12 CO are < 1/25 and < 1/100 respec- 
tively. This hot component is needed to fit correctly the 



high-J CO lines, as the single 400 K component of Cox et 
al. (1996) does not seem to be sufficient in our view. 

(ii) As this fast wind slows down when shocking with the 
outer parts of the envelope, we have considered intermedi- 
ate layers of larger sizes. Parameters are from Cox et al. 
(1996, 2000): a layer of 500 K, 8=2" , 100 kms" 1 expan- 
sion velocity, and density of 5 10 7 cm -3 . The 12 CO column 
density has decreased to 7 10 17 cm~ 2 and the HCN/CO 
abundance ratio is < 1/30. An additional component is 
required where the wind has slowed down to 60 kms -1 , 
with a temperature of 400 K, 9—3", density of 5 10 7 cm~ 3 
and 7V( 12 CO )= 5 10 17 cm" 2 . 13 CO and HCN have the 
same relative abundance ratios than in the previous layer. 

(iii) A moderate velocity wind region (Sahai et al. 1998) 
with 8—15", temperature of 250 K, density of 5 10 7 cm~ 3 , 
velocity of 40 kms" 1 , and 12 CO column density of 4.5 10 16 
cm~ 2 . The relative abundances for the other species are 
the same as for previous layers. 

(iv) The AGB remnant CSE where low-excitation CO 
emission extends beyond 40" (Cox et al. 1996, Yamamura 
et al. 1996) and that we have modelled as a layer of 25" at 
50 K, expansion velocity of 20 kms" 1 , iV( 12 CO )= 3.5 10 16 
cm . 

The global fit is shown in Figure 3 and reproduces rea- 
sonably well the LWS-ISO grating observations of CRL 
2688. We also applied the same model to the mm and 
submm CO lines. Components (iii) and (iv) mainly fit 
the low-J CO lines (millimeter and submillimeter lines). 
Observed and predicted intensities are in good agreement 
(see Table 7). 

5.3. Changes since the previous stage 

With respect to the previous AGB stage, represented 
by IRC+10216, we see that the molecules emitting in the 
far-IR are basically the same, but with different relative 
intensities. This is probably due to the fact that the in- 
nermost region of the AGB remnant envelope has been 
already pushed away in these initial stages of the evolu- 
tion towards the PPN phase. Hence, all the ro-vibrational 
lines of HCN have disappeared, not necessarily due to a 
depletion of this molecule but to the lack of sufficiently 
high temperatures and densities. This fact, together with 
the lack of [OI] and [CiT], makes CO the main coolant of 
the envelope. The observational evidence of winds faster 
than those in the previous stage indicates that the gas in 
the inner layers is now heated by shocks, where the high-J 
CO lines excitation is favored. The presence of polyynes, 
as shown by the ISO mid-IR spectrum of this object (Cer- 
nicharo et al., 2001a), indicates that the envelope has also 
started an important chemical evolution when compared 
to IRC+10216. The high velocity winds and the increas- 
ing UV field from the central object are at the origin of 
this new chemistry, and this will be enhanced in the next 
stage (CRL 618). 

6. A PPN OBJECT: CRL 618 

CRL 618 is one of the few clear examples of an AGB 
star in the transition phase to the Planetary Nebula stage. 
There is evidence that this object must be a <~ 200 years 
old PPN (Bujarrabal et al. 1988). It has a compact HII 
region created by a hot central star (30 000 K, spectral 
type B0, Justtanont et al. 1997). CRL 618 is seen as a 
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bipolar nebula at optical, radio and infrared wavelengths 
(Carsenty & Solf 1982, Bujarrabal et al. 1988, Cernicharo 
et al. 1989, Neri et al. 1992, and Hora et al. 1996). The ex- 
pansion velocity of the outer envelope is around 20 kms -1 . 
CO observations by Cernicharo et al. (1989) evidenced the 
presence of a high-velocity outflow (^200 kms -1 ). The 
UV photons from the central star and this high velocity 
wind perturb the circumstcllar envelope (CSE), produc- 
ing shocks and photodissociation regions which modify the 
physical and chemical conditions of the gas (Herpin & Cer- 
nicharo 2000; Cernicharo et al. 2001a&b). 

The molecular content of CRL 618 (Cernicharo et 
al. 2002, in preparation) is characterized by a large 
abundance of cyanopolyynes, acetylenic chains, methane, 
mcthylpolyynes and benzene (Cernicharo et al. 2001a&b): 
this PPN is a very eficient organic chemistry factory. The 
detection of long carbon chains in both CRL 618 and CRL 
2688 suggests that the input of a UV radiation field from 
the evolving central object or/and the action of high veloc- 
ity winds on the slowly expanding AGB remnant modify 
the CSE chemistry allowing the formation of new and more 
complex molecules. 




187 188 189 190 



Fig. 4 -The HCN J=18-17 line emission detected by the LWS 
grating in CRL 2688 at a resolution of 0.6 /im (S/ ~ 4). 

6.1. Low- excitation CO lines 

Two outflows can be seen in the CO line profiles of CRL 
618 (see Table 2). The main outflow is centered at —22.0 
kms -1 with an expansion velocity of roughly 20 kms -1 . 
The optically thick CO J=l— line clearly absorbs the 
continuum from —60 to —40 kms -1 coinciding with the 
main outflow terminal velocity. Other lines (higher tran- 
sitions) show an absorption dip (self-absorption). In CRL 
618, many molecular lines show P-Cygni profiles at 3 and 
2 mm (Cernicharo et al. 2002). These lines, which corre- 
spond to very high energy ro-vibrational levels of HC3N 
and HC5N are formed in the PDR of the high density torus 
surrounding the central object. 

The HVW (around 200 kms -1 ), associated with a bi- 
conical outflow, is also observed in many other molecular 
transitions (see also Cernicharo et al. 1989, Cernicharo et 
al. 2002) but with a smaller velocity span. The CO wings 



corresponding to this outflow have a smooth triangular 
shape and are present in all the millimeter and submil- 
limctcr lines. Their intensity increases with J, Rj going 
from 0.09 for J=l — to an estimated value of more than 
2.6 in the J =7— 6 line. Therefore, the emission in the high 
J lines (next paragraph) is basically dominated by this 
high velocity wind component. 

6.2. Far-IR emission 

The rotational lines of 12 CO (J=14-13 to J=41-40), 
13 CO (J=14-13 to J=19-18), HCN, HNC, [OI] (at 
63.170 and 145.526 /um) and [CI I] ( 2 p 3/2 - 2 P1/2 tran- 
sition at 157.74 /j,m) are observed in CRL 618 (see Fig. 
3). Several lines of OH and H2O are also detected. How- 
ever, no H2O emission has been detected at 183 GHz. As 
the 22, 183 and 321 GHz water mascr lines are collision- 
ally excited (radiative effects are small) in warm dense gas 
(Neufeld & Melnick 1990, 1991), this non-detection may 
indicate that water is not in the wind region, but rather 
in the torus. 

The model used for analyzing the LWS observations of 
CRL 618 was presented in Herpin & Cernicharo (2000). 
Briefly, it consists of a central high density torus with 
its inner part being a Photo-Dissociated Region, an ex- 
tended AGB remnant envelope and high velocity bipolar 
outflows. The physical parameters of the different layers 
are given in Table 5. The AGB remnant parameters are 
derived from Yamamura et al. (1994) and Hajian et al. 
(1996). The HVW characteristics are from Cernicharo et 
al. (1989), Neri et al. (1992) and Martin-Pintado et al. 
(1995). Justtanont et al. (2000) fit the ISO CO lines with 
a 700 K component (n(H2)=5 10 7 cm -3 ). Temperatures 
of the torus layers have been varied in order to get the 
best fit of the data. Our model, once applied to the low- 
excitation CO lines, disagrees with the observed J=3— 2 
to 7—6 intensities. This may be mainly explained by the 
completely non-spherical geometry of the source, and by 
the complex velocity field across the envelope. Moreover, 
as the CO molecules are reprocessed, excitation is more 
complex now and cannot be reproduced by our relatively 
simple model. 

6.3. Changes since the previous stage 

As in the previous very young PPN phase, shock chem- 
istry plays an important role, as even faster outflows are 
detected. The now higher effective temperature of the 
central star induces a UV-based chemistry dominated by 
the photolysis of key molecules and the molecular con- 
tent has started to be reprocessed. This fact is proved by 
the observed changes in molecular abundances (cf. HCN, 
HNC), the huge abundance of cyanopolyynes, the pres- 
ence of polyynes and benzene (Cernicharo et al., 2001a&b) 
and of O-bearing molecules (Herpin and Cernicharo 2000). 
Ionic ([Ci7]) and atomic fine structure lines ([O/]) are also 
detected as a result of the photodissociation of CO. At this 
stage, the main cooling of the gas occurs via CO high-J, 
[CI I] and [OI] lines. 

7. A YOUNG PN OBJECT: NGC 7027 

NGC 7027 is the most evolved object in this study. It 
is a very young PN, having left the AGB only 10 3 years 
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ago (Volk & Kwok 1997). The estimated temperature of 
the central star is around 200 000 K (Latter et al. 2000). 
The gas in the inner 10" radius of the envelope has been 
completely ionized and a strong continuum emission arises 
from it. This region is revealed in the far-infrared by 
fine-structure lines (Liu et al. 1996). The PDR at the in- 
terface between the ionized and neutral regions is mainly 
traced by the pure rotational lines of CH + (Cernicharo et 
al. 1997). The H 2 emission at this interface is found to 
be consistent with excitation by UV photons (Cox et al. 
1997). 

7.1. Low-excitation CO lines 

All the CO lines are asymmetrical (see Fig.l), reflecting 
the complex kincmatical structure of this PN. The main 
CO line emission is centered at Vlsr ^27 — 30 kms" 1 , 
with an expansion velocity of 15 kms" 1 for the J=l— 0, 
2—1, 3—2 and 4—3 lines, and 20 kms" 1 for the 6—5 and 
7—6 lines (see Table 3). The higher velocity extent of 
the last two transitions probably indicates the presence of 
a faster wind at smaller radii (3" — 5"). Jaminet et al. 
(1991) found 

v e xp ~ 15 kms 1 (with an angular exten- 
sion of 5" — 15") for the 2—1 and 3—2 lines in agreement 
with our data for those transitions. They also have de- 
rived a velocity dispersion as high as 3 kms -1 at small 
radii (<5"), and more than 1 kms -1 in the outer parts 
(>5"). An absorption feature at vlsr ~ 10 kms" 1 , due 
to self- absorption (Jaminet et al. 1991), is detected in the 
low-J lines (sec Fig. 1). 

Although the J=l— line of CO is certainly slightly con- 
taminated by the 38a recombination line of atomic hy- 
drogen (shifted only by —8.3 kms -1 with respect to its 
center), the line profiles show also some emission at high 
velocity (mainly in the red part) that could be the signa- 
ture of a faster outflow, with an estimated velocity extent 
of 45 kms" 1 . This outflow component was not reported 
by Jaminet et al. (1991). They only found the extended 
outflow with v exp ~ 15 kms -1 and a wind at a smaller ra- 
dius with v exp ~ 23 km s" 1 , which may be a flow driven by 
the pressure of the ionized region, and which is probably 
the place where the 6—5 and 7—6 lines are mainly excited 
(see v m ain in Table 3). The flow at higher velocity (^47 
kms -1 ) may be the relic of a High Velocity Wind ejected 
during the PPN phase. Nevertheless, the importance of 
the wing emission is quite low as the ratio Ri is less than 
0.1. 

7.2. Far-IR emission 

The J=14-13 to J=23-22 12 CO lines are detected 
in the LWS grating spectrum of NGC 7027 and some 
13 CO lines are also marginally detected. Despite its car- 
bon rich chemistry, several lines of OH are clearly identi- 
fied (Liu et al. 1997), while no H2O emission is observed 
(Cernicharo et al. 1997). In addition to the molecu- 
lar emission, strong atomic and ionic features typical of 
a PN spectrum are also seen. They include [OI], [CI I], 
[Nil], [Nil I] and [0777"] lines. The pure rotational lines of 
CH + discovered by Cernicharo et al. (1997) are prominent 
through the far-infrared spectrum. Features at 149.18 and 
180.7 /im correspond to rotational lines of CH (see Liu et 
al. 1997). A new feature of this species at 115.9 [im is ten- 
tatively detected. However, no HCN or HNC emission is 



seen, which supports the small HCN abundance (9 10~ 9 ) 
derived by Deguchi et al. (1990). The main cooling agents 
arc fine-structure atomic and ionic lines. 

To fit the far-infrared spectrum of NGC7027 we have 
considered 3 components (see Table 6): (i) the atomic re- 
gion; (ii) a high density layer; (iii) the AGB remnant. Sizes 
of the different regions are derived from Yan et al. (1999) 
and Jaminet et al. (1991). 

The atomic region (i) has been subdivided into two lay- 
ers of 1000 K and 800 K and 6= 5" and 6" respectively 
(Cernicharo et al. 1997, Liu et al. 1996, Hasegawa, Volk & 
Kwok 2000). The density is 5 10 7 cm" 3 , and the expan- 
sion velocity 15 kms" 1 . In the hottest layer, only [07] and 
CH + are present (for now we have not considered [C77] 
and [0777], in the model although they are also present). 
In the external layer, we introduce 12 CO , 13 CO and OH. 
The considerable abundance of OH suggests that H 2 was 
probably present in the previous stage of evolution (see 
CRL 618), and, even if not detected, may be still present 
here, but with a low abundance. Thus we placed H 2 in 
region (i) with an upper limit of abundance relative to H 2 
of 1.5 10" 7 , according to the present data. This upper 
limit is 50% lower than the one obtained by Volk & Kwok 
(1997). The best fit to the data gives the following abun- 
dances (relative to 12 CO ) in this layer: < 1/40, 1/20, 1/80 
and 875 for 13 CO , OH, CH+ and [07] respectively. Com- 
paring with analyses performed by other authors we point 
out that the H 2 densities used in our model are higher 
than those of Liu et al. (1996), who took values between 
10 5 and 10 6 cm -3 for the high-J CO lines coming from the 
hot region (1000—700 K). Also, the velocity of the expand- 
ing ionized region derived by Jaminet et al. (1991) is 17.6 
kms" 1 , while the estimated CO velocity wind by Sopka 
et al. (1989) is 17 kms" 1 . Note that the two components 
of the atomic region produce similar lines, except that the 
hottest one only contains [07] and CH + . 

In the high density layer (ii) we consider a temperature 
of 270 K, 6— 9", expansion velocity of 15 kms" 1 and a 
density of 10 6 cm" 3 (Yan et al. 1999, Justtanont et al. 
2000). Volk & Kwok (1997) applied for this layer very 
similar values (n H2 =9 10 5 cm" 3 and 7 = 220 K). The 
fit results in a 12 CO column density of 2.1 10 17 cm" 2 and 
[ 12 CO/ 13 CO]<30. 

The AGB remnant region (iii) has an extension go- 
ing from 10" to 20", the gas temperature decreases from 
100 K to 50 K, and the density from 3 10 5 to 10 4 cm" 3 
(Hasegawa, Volk & Kwok 2000). The expansion velocity is 
15 kms" 1 (see Sect. 6.1 and Volk & Kwok 1997, Jaminet 
et al. 1991). The colder layer (50 K) is important here 
only to explain the low-J CO emission. The 12 CO column 
density is 1.5 to 1.7 10 17 cm" 2 and the [ 12 CO / 13 CO ] ra- 
tio is always < 30. Note that CO emission extends over a 
region larger than 40"(Masson et al. 1985). 

Our lower limit of the 12 C/ 13 C ratio (30) is compatible 
with the value derived by Kahane et al. (1992) for this 
object (65) using much more sensitive millimeter observa- 
tions of several carbon-bearing molecules. The whole CO 
column density derived by Thronson (1983) is larger than 
2 10 17 cm" 2 , in good agreement with our values. 

The results of our model when applied to the low- 
excitation CO lines are satisfactory for J=l— to J=4— 3, 
but the predicted intensities are too high for the J=6— 5 
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and 7—6 lines. The reason might be the complex kinemat- 
ical structure of this PN. 

7.3. Changes since the previous stage 

The main characteristic of this young PN compared with 
previous stages is the almost disappearance of fast molec- 
ular winds. The hot central object induces a strong UV 
field leading to a UV-dominated photochemistry in the 
envelope. As a consequence, molecular abundances have 
been strongly modified since the PPN phase. Water va- 
por and HCN have almost disappeared, while other species 
such as CH + have been efficiently formed. The cooling ef- 
fect of CO lines is much weaker and now the gas is mainly 
cooling via ionic and atomic fine structure lines. 

8. DISCUSSION 

• It is assumed that CRL 2688 has just left the AGB 
phase: the central object is not yet sufficiently hot to pro- 
duce a significant ionization or photodissociation of the 
gas. Hence, the far-IR molecular emission ( 12 CO , 13 CO , 
HCN) is well reproduced by our model. Nevertheless, the 
detection of the polyynes C4H2 and C6H2 in this object in- 
dicates that some chemical processing has already started, 
an important difference as respects to IRC+10216 where 
these species are not detected. 

In addition to the low expanding AGB CSE remnant, 
this PPN exhibits an inner and faster outflow which runs 
into that AGB envelope. This is one of the main charac- 
teristics of an object beginning its transition to the PN 
stage. At this point, emission from the molecular gas is 
the main cooling process. 

The 13 CO emission is not detectable at the low sensitiv- 
ity of the LWS grating spectrum. For a 12 CO abundance 
of 6 10~ 4 the upper limit of 13 CO abundance would be 
< 2.2 10~ 5 . This value is half of what was found by Ya- 
mamura et al. 1996 (4 10~ 5 - 7.5 10~ 5 ), but is similar 
to the 13 CO abundance found in an AGB star such as 
IRC+10216 (1-3 10~ 5 , Cernicharo et al. 1996). However, 
Kahane et al. (1992) derived more accurate isotopic ra- 
tios, 3.7 10~ 5 and 7.5 10" 5 , respectively for IRC+10216 
and CRL 2688, showing an increase by a factor of 2. For 
Jaminet et al. (1992), CNO nuclear processing in the star 
injects material into the stellar wind: in the slow wind, 
[ 12 CO / 13 CO ] = 20, while in the fast wind they find a 
ratio of 5. Because of these values, these authors placed 
CRL 2688 among the J-type carbon stars, which have a 
considerably lower 12 CO / 13 CO ratio than carbon stars 
(40—80). However, we prefer to argue that CRL 2688 has 
just begun to evolve to the PN stage; the 12 CO material 
has just began to be reprocessed, as proved by the appar- 
ent increase of the 13 CO abundance between IRC+10216 
and CRL 2688. 

The HCN/CO ratio in IRC+10216 is 0.1, much larger 
than in the high velocity winds of the post-AGB en- 
velopes discussed here. According to our model, we find a 
lower HCN abundance (relative to 12 CO ) in the fast wind 
(< 6 10~ 6 ) than in the slow wind (< 2 10~ 5 , i.e., similar 
to IRC+10216). Sopka et al. (1989) found 3.3 10~ 6 in the 
fast wind while Jaminet et al. (1992) found 7 10" 7 . 

HNC is not detected in this object. This molecule is 
mainly formed by ion-molecule reactions, and in this post- 
AGB envelope, physical conditions are therefore not very 



suitable for a considerable HNC production. Moreover, 
Yamamura et al. (1996) and Sopka et al. (1989) derived 
[HCN/HNC] ~ 21 and 160 respectively. According to this 
fact, the HNC emission is largely under our detection limit. 

• In CRL 618, Herpin & Cernicharo (2000) have shown 
that O-bearing species, H 2 and OH, are produced in the 
innermost region of the circumstellar envelope. Also Cer- 
nicharo et al. (2001 a&b) have detected in this object the 
poli-acetylenic chains C4H2 and C6H2, methyl-polyynes, 
and benzene. The UV photons from the central star pho- 
todissociate most of the molecular species produced in the 
AGB phase and allow a chemistry dominated by stan- 
dard ion-neutral reactions. O-bearing species are formed, 
and also abundances of C-rich molecules such as HCN and 
HNC are modified (see Table 7) . At the high temperatures 
of the PDR, OH can be formed by endothermic reactions 
between O and H 2 , with O released by CO photodissoci- 
ation, and then can produce H 2 through reactions with 
H 2 . H 2 and OH can also be produced from the dis- 
sociative recombination of H 3 + if this ion is abundant 
enough. 

• NGC 7027 has a very hot central star. Because of the 
high UV flux (10 4 < G < 10 5 , Burton et al. 1990), there 
is a strong presence of atomic lines. UV photons from the 
central star produce a PDR and the gas cools mainly via 
fine-structure atomic lines. In the PDR, the total number 
of CO molecules is around 2% of that of ionized carbon, i.e. 
most of the CO molecules have been photodissociated (Liu 
et al. 1996). As the ionized region is not modeled in this 
work, our results concerning the atomic region have to be 
carefully taken: in the high electron density environment 
of NGC 7027, ionic fine-structure lines are strongly sup- 
pressed by collisional de-excitation (Liu et al. 1996, Keyes 
et al. 1990). Electron impact may be important in excit- 
ing the CH + lines and would lower the CH + abundance 
(Cernicharo et al. 1997). In the 1000 K layer, the forma- 
tion of CH+ may be achieved via C+ + H 2 ->■ CH+ + H 
(activation energy of ~ 4000 K). The formation rate of 
CH + (Sternberg & Dalgarno 1995) could be high enough 
to produce column densities similar to those derived in this 
work. Due to the high temperatures, densities and UV ra- 
diation field in the PDR, the formation of CH + leads to 
the creation of CHj and CH^, whose dissociative recom- 
bination will form CH very efficiently at a temperature of 
800 K (Sternberg & Dalgarno 1995); Liu et al. (1997) de- 
rived [CH/CH+]- 0.2. The observed CH lines are excited 
by collisions with atomic and molecular hydrogen. 

The presence of [OI] and CH+ in the 1000 K layer of 
the atomic region, and the low abundance of CO in the 
800 K layer, indicate that the CO molecules have been 
largely reprocessed there through UV photons. Note that 
the chemical model of Hasegawa, Volk & Kwok (2000) 
shows that the CO observed in these objects consists of 
new formed molecules. Nevertheless, the AGB circumstel- 
lar gas has not been totally reprocessed: the ionized and 
photodissociation regions still constitute a relatively small 
but significant fraction of the total mass of the circumstel- 
lar material around NGC 7027. 

The study of CRL 618 showed that H2O and OH can be 
efficiently produced. However, these O-bearing species are 
also quickly reprocessed. This is seen in the next stage, 
represented by NGC 7027. The equilibrium models (Mc- 
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Cabe, Connon, Smith & Clegg 1979) predict essentially 
no water in the carbon-rich conditions of the NGC 7027 
envelope. Volk & Kwok (1997) argue that even under non- 
equilibrium chemical conditions, only a very small fraction 
of oxygen will be converted to water here. The non detec- 
tion of H 2 and the detection of OH in NGC 7027 in- 
dicate that H 2 molecules formed in the previous stage 
may have been almost entirely reprocessed. In the NGC 
7027 environment, the simplest explanation for the water 
vapor disappearance is its transformation into OH and H 
by UV photo-dissociation. This will explain the increase 
of OH abundance by a factor 60 since the previous PPN 
stage represented by CRL 618. Concerning HCN, prob- 
ably most molecules have been photodissociated into H 
and CN (CN was detected by Bachiller et al. 1997 with 
a strong abundance, CN/HCN=9). Photo-dissociation of 
CO, HCN, and other carbon-bearing molecules contribute 
to the increase of carbon ions abundance. 

9. CONCLUSION 

The study of the selected 3 objects in rapid evolution 
from the AGB to the PN stage clearly shows the cru- 
cial importance that (i) the interaction between fast and 
slow winds; and (ii) the increase of UV flux as the central 
object evolves toward the white dwarf state have on the 
chemistry. The spectroscopic evidence of the increasing 
UV flux comes from the appearance of atomic and ionic 
lines in CRL 618, which later dominate the far-IR spec- 
trum of NGC 7027. Atomic and ionic lines tend to appear 
(at ISO's sensitivity) when the central star is hotter than 
10000 K (Fong et al. 2001; Castro-Carrizo et al. 2001). 
On the other hand, the spectroscopic signature showing 
that shocks become less important as the evolution goes 
on, comes from the wind velocity decrease seen in the CO 
profiles (this must be confirmed by further studies on more 
objects and only refers to the neutral gas). The strongest 
shocks occur just after leaving the AGB when the central 



star is ejecting large amounts of material in a very fast 
wind (the case of CRL 2688). The AGB remnant enve- 
lope is being shocked by an inner, faster wind developed 
in the PPN stage (i.e., the 200 kms" 1 wind in CRL 618). 
Most of the 12 CO , 13 CO and HCN emission is produced 
in these shocks. The fast increase of the stellar temper- 
ature, will produce a new UV-dominated chemistry when 
T e // >30000 K. These new conditions (UV photons and 
shocks) will deeply modify the constitution of the inner 
parts of the envelope. Indeed, in CRL 618, O-bearing 
molecules (H 2 and OH, Herpin and Cernicharo 2000) 
and relatively complex organic molecules (Cernicharo et 
al., 2001a&b) appear. Furthermore, most of the CO and 
HCN will be entirely reprocessed in the PDR leading to 
strong HNC emission. At this point, CO and [OI] atomic 
lines are the dominant coolants. 

As the star reaches the PN stage, the strong fast molec- 
ular winds have disappeared, and slow expanding layers 
constitute the PN envelope around a large and hot atomic 
region. Most of the old AGB material has been repro- 
cessed. The spectrum is now dominated by atomic and 
ionic lines. New species such as CH + and CH appear. 
There is only weak HCN emission, as the molecules may 
have been broken into H and CN. More interesting is the 
disappearance of H 2 , which has probably also been re- 
processed, and is only a relatively abundant molecule in 
the intermediate C-rich PPN stage. 
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TABLE 1 

Tabic of the parameters for the different flows in CRL 2688. The uncertainty on v main and v W i ng is 0.3 krns" 1 . Ri is the 
wing and main integrated intensities ratio. 









CRL 2688 






CO parameters 


1-0 


2-1 


3-2 


4-3 


6-5 


7-6 


Vmain (kins' 1 ) 


19.3 


20.0 


23.9 


24.1 


23.4 


23.3 


Vwing (km S ) 


50.0 


52.0 


50.7 


40.3 


41.6 


33.3 


Ri 


0.08 


0.08 


0.11 


0.16 


0.24 


0.24 


ARi 


(0.01) 


(0.01) 


(0.01) 


(0.02) 


(0.02) 


(0.03) 



TABLE 2 

Table of the parameters for the different flows in CRL 618. The uncertainty on v ma in and v wing is 0.3 krns" 1 . Ri is the 
wing and main integrated intensities ratio. For the 6—5 and 7—6 lines, the parameters Vwing and Ri are lower limits due 
to the unsufhcicnt bandwidth available for these observations. 









CRL 618 




CO parameters 


1-0 


2-1 


3-2 


6-5 


7-6 


Vmain (klTlS' 1 ) 


19.7 


19.6 


18.9 


17.1 


18.8 


Vwing (kms _1 ) 


202.5 


204.5 


178.5 


>76.0 


> 195.0 


Ri 


0.09 


0.22 


0.45 


>0.72 


>2.6 


ARi 


(0.01) 


(0.02) 


(0.03) 


(0.05) 


(0.3) 



TABLE 3 

Table of the parameters for the different flows in NGC 7027. The uncertainty on v ma i n and v W i ng is 0.3 kms' 1 . Rj is the 
wing and main integrated intensities ratio. 









NGC 


7027 






CO parameters 


1-0 


2-1 


3-2 


4-3 


6-5 


7-6 


Vmain (kms' 1 ) 


14.6 


14.3 


16.0 


18.3 


20.0 


20.0 


V w ing (kms _1 ) 


47.2 


43.2 


30.0 


25.1 


24.7 


23.8 


Ri 


0.04 


0.09 


0.07 


0.07 


0.06 


0.09 


AR! 


(0.01) 


(0.01) 


(0.02) 


(0.02) 


(0.02) 


(0.02) 



TABLE 4 

Table of the column densities (in cm -2 ) obtained from our model for all molecules present in each region of CRL 2688 (6 
defines the solid angle of each considered layer). 



Molecules 


"CO 


3 CO 


HCN 


FAST WIND FROM THE STAR (i) 

T=800 K, 9 = 0.3", v = 200 kms' 1 , n(H 2 ) =5 10 7 cm' 3 
T=500 K, 6 = 2", v = 100 kms' 1 , n(H 2 ) =5 10 7 cm' 3 
T=400 K, e = 3", v = 60 kms" 1 , n(H 2 ) =5 10 7 cm" 3 


2 10 18 
7 10 17 
5 10 17 


<8 10 16 
<2.8 10 16 
<2 10 16 


<2 10 16 
<2.3 10 16 
<1.7 10 16 


MEDIUM WIND FROM THE STAR (ii) 

T=250 K, 9 = 15", v = 40 kms" 1 , n{H 2 ) =5 10 7 cm' 3 


4.5 10 16 


<1.8 10 15 


<1.5 10 15 


AGB remnant (iii) 

T=50 K,e = 25", v = 20 kms -1 , n(H 2 ) =5 10 5 cm" 3 


3.5 10 16 


<1.4 10 15 


<1.2 10 15 
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TABLE 5 

Tabic of the column densities (in cm -2 ) obtained from our model for all molecules present in each region of CRL 618 (8 
defines the solid angle of each considered layer). 



Molecules 


"CO 


13 co 


HCN 


HNC 


o-H 2 


OH 


[OI] 


TORUS (including PDR) 
T=250 K, = 1.5" 


10 19 


5 10 17 


10 16 


10 16 


3 10 17 


8 10 15 




v = 20 kms" 1 , n{H 2 ) =5 10 7 cm" 3 
















T=800 K, = 1.1" 


6 10 17 


3 10 16 












v = 20 kms" 1 , n(H 2 ) =5 10 7 cm" 3 
















T=1000 K, = 0.6" 


10 19 


5 10 17 










4.5 10 19 


v = 20 kms" 1 , n(H 2 ) =7 10 7 cm" 3 
















HVW 
















T=200 K, = 1.7" 


5 10 18 


2.5 10 17 


3 10 17 


2 10 17 








v = 50 kms" 1 , n(H 2 ) = 10 7 cm" 3 
















T=1000 K, 9 = 1.5" 


2 10 16 


1 10 15 


2 10 15 


10 14 








v = 200 kms" 1 , n(H 2 ) = 10 7 cm" 3 
















AGB remnant 
















T=50 K,0 = 10" 


7 10 17 


3.5 10 16 


7 10 16 


7 10 14 








v = 20 kms" 1 , n{H 2 ) = 10 5 cm" 3 

















TABLE 6 

Table of the column densities (in cm" 2 ) obtained from our model for all molecules present in each region of NGC 7027 (6 
defines the solid angle of each considered layer). 



Molecules 


"CO 


13 CO 


o-H 2 


OH 


[OI] 


CH+ 


ATOMIC REGION (i) 

T=1000 K, = 5" 

v=15 kms" 1 , n{H 2 ) =5 10 7 cm" 3 

T=800 K, = 6" 

v = 15 kms" 1 , n(H 2 ) =5 10 7 cm" 3 


8 10 15 


2.7 10 14 


<2.6 10 13 


4 10 14 


2 10 19 
7 10 18 


9 10 13 
10 14 


HIGH DENSITY LAYER(ii) 

T=270 K, = 9" 

v = 10 kms" 1 , n(H 2 ) = 10 6 cm" 3 


2.1 10 17 


7 10 15 










WIND REGION (iii) 

T=100 K, 9 = 10" 

v = 15 kms" 1 , n(H 2 ) =3 10 5 cm" 3 

T=50 K, = 20" 

v = 15 kms" 1 , n(# 2 ) = 10 4 cm" 3 


1.5 10 17 
1.7 10 17 


5 10 15 
5.7 10 15 











TABLE 7 



Observed and predicted intensities for the mm and submm CO lines for each source. 



Source 


CO lines 




1-0 


2-1 


3-2 


4-3 


6-5 


7-6 


CRL 2688 


WJy) 


30 


70 


271 


373 


640 


760 




Rw(Jy) 


10 


30 


296 


429 


650 


790 


CRL 618 


io6 3 ( J y) 


23 


47 


79 




320 


232 




Imod(Jy) 


21 


66 


197 


400 


742 


994 


NGC7027 


Lb S (Jy) 


36 


70 


364 


516 


520 


840 




Imod(Jy) 


24 


83 


475 


498 


893 


1061 



TABLE 8 

Table of the molecular, atomic and ionic abundances (relative to 12 CO ) in CRL 2688, CRL 618 and NGC 7027 as seen 
by the model. 



Species 


CRL 2688 


CRL 618 


NGC7027 


13 CO 


< 1/25 


1/20 


£1/30 


HCN 


< 1/30 - 1/100 


1/10-1/1000 




HNC 




1/10-1/1000 




H 2 




1/25 


< 1/650 


OH 




1/1250 


1/20 


[OI] 




4.5 


875 


CH+ 






1/80 



